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Abstract 

Distribution of the conductance P(g) at the critical 
point of the metal-insulator transition is presented for 
three and four dimensional orthogonal systems. The 
form of the distribution is discussed. Dimension de- 
pendence of P(g) is proven. The limiting cases g — > oo 
and 5^0 are discussed in detail and relation P(g) — > 
in the limit g —> is proven 

PACS numbers: 71.30.+h, 71.23.-k, 72.15.Rn 

As the conductance g in disordered systems is not 
the self-averaged quantity, the knowledge of its prob- 
ability distribution is extremely important for our un- 
derstanding of transport. This problem is of special 
importance at the critical point of the metal-insulator 
transition [|l]. While the distribution of the conduc- 
tance in the metallic phase is known to be Gaussian in 
agreement with the random-matrix theory || and the 
localized regime is characterized by the log-normal dis- 
tribution of g the form of the critical distribution 
remains still unknown. Among the problems which are 
not solved yet we mention e.g. the number of param- 
eters which characterize distribution, the existence of 
huge fluctuations of the conductance, and the form of 
P(g) for small values of g. 

Several attempts has been made to characterize con- 
ductance distribution at the critical point. Using the 
Migdal - Kadanoff renormalization treatment, huge 
conductance fluctuations has been predicted in ||] . The 
same conclusion was found also in systems of dimen- 
sion d = 2 + e. In the limit e « 1 the form of the 
distribution P{g) was found analytically [0. However, 
numerical studies of disordered 3D system || indicated 
that it is not possible to generalize these analytical con- 
clusions for realistic 3D systems (e = 1). 

The form of P(g) for 2D symplectic models was 
found in || [7]. Recently, P(g) has been studied also 
for system in magnetic field, both in 3D B and in 2D 



|9|. The main conclusion of these studies is that the 
symmetry of the system influences the form of the dis- 
tribution at the critical point more strongly than in 
the metallic or localized regime. Nevertheless, P(g) is 
invariant with respect to the choice of the microscopic 
model within the same universality class . 

Studies of the statistics of the conductance have 
their counterpart in the analysis of the level statistics 
s = JSj+i — Ei of the eigenvalues of Hamiltonian [fTT|| . 
The critical distribution P(s) is also the subject of in- 
tensive studies within last year s Il2[ . In particular, its 
dependence on the symmetry [|l3|, and dimension fl4f| 
have been studied numerically. 

In this Letter we present new numerical data for the 
3D and 4D Anderson model (orthogonal ensembles). 
Although data prove the dimension dependence of the 
distribution, their enables us to discuss the common 
features of the critical distribution. In particular, we 
prove that P(g) decreases more quickly than exponen- 
tially for large g. This assures that there are no huge 
fluctuations of the conductance, discussed in ||. We 
prove also that P(g) — ► in the limit of g — * 0. 

We calculated the conductance as 



!J 



= Tr th = ^cosh~ 2 (>j/2) 



(1) 



where quantities Z{ determine eigenvalues of the trans- 
mission matrix t't. Details of the method have been 
published elsewhere ||. For a given system size L, the 
probability distribution of g has been calculated from 
an ensemble of iV s t a t samples. The list of used ensem- 
bles together with mean and variances of g are given 
in Table 1. 

The last column of Table 1 presents parameter 
(zi), which corresponds to the parameter A introduced 
in the finite size scaling theory by MacKinnon and 
Kramer as {zi) — in the quasi-one dimen- 
sional limit L d ^ 1 x Lt, L t >> L. When neglecting the 
smallest system size, our data confirm the i-invariance 
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Figure 1: Probability distribution of log g /(g) (left) 
and g/(g) (right) for 4D (full symbols) and 3D (open 
symbols) Anderson model. For comparison, we plot 
also data for 2D (symplectic) Ando model. The last 
exhibit the best convergence for both small and large 
values of g. For meaning of symbols, see Table 1. Solid 
line is Poisson distribution P(g) — exp —g/(g). 

of (z\) as well as of (g) and (log g) and their standard 
deviations. Owing to higher critical disorder, (z\) is 
larger in 4D than in 3D. This guarantees that the finite- 
size effects disappear more quickly in 4D. Therefore, in 
spite of the fact that computer facilities limited the 
system size to L < 8 for d — 4, obtained data provide 
us with the relevant information about all parameters 
of interest. 

We presented in Table 1 both mean values of g and 
log g to underline the common features of 3D and 4D 
distribution: the variance of log g is of order of its mean 
value. This relation is typical for localized state. On 
the other hand, standard deviation of g is also ~ (g). 
Its value for 3D samples, 0.334, is smaller than the 
same quantity calculated for 3D in the metallic regime 
1- 

Numerical data for P(g) are presented in Figure 
1. They confirm that the critical distribution of g is 
system-size independent, in agreement with previous 
studies. Fig. 1. shows also that P(g) depends on the 
dimension of the system within the same symmetry 
class. Although the distribution has the same form for 
3D and 4D ensembles, it becomes broader for higher d: 
the probability to find g « (g) or g » (g) growths 
with dimension. This is due to higher critical disorder, 
which causes that electronic state posses more features 



Table 1: Review of ensembles studied in the present 
work. L-size of the d-dimensional cube, N a t a t. number 
of samples in a given ensemble, var g = (g) 2 — (<? 2 ), (z\) 
is mean of the smallest of z's. Data for 3D AM are in 
good agreement with |)| (up to the spin degeneracy 
factor 2). 
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3D Anderson model: W c 


« 16.5 






6 


o 


20.000 


0.375 


0.324 


-1.481 


1.344 


2.901 


8 


< 


20.000 


0.400 


0.333 


-1.384 


1.251 


2.803 


10 


> 


10.000 


0.410 


0.337 


-1.347 


1.229 


2.770 


12 





5.000 


0.421 


0.340 


-1.302 


1.199 


2.724 


14 


A 


2.500 


0.416 


0.338 


-1.306 


1.122 


2.725 


18 


V 


500 


0.418 


0.329 


-1.279 


1.083 


2.717 






4D Anderson model: W c « 


34.5 






4 




22.000 


0.190 


0.247 


-2.569 


2.301 


4.130 


5 


V 


30.000 


0.229 


0.270 


-2.275 


2.006 


3.838 


6 


□ 


15.000 


0.225 


0.269 


-2.291 


2.054 


3.852 


7 


A 


7.000 


0.239 


0.275 


-2.193 


1.971 


3.748 


8 




200 


0.227 


0.274 


-2.188 


1.692 


3.790 



of the localized state than that of the metallic one (re- 
maining critical). This is in agreement with studies of 
the level statistics in 4D 

The small- g behavior of P(g) can be estimated from 
Fig 1. Instead of P(g/(g)), we plot in the left side of 
Fig. 1. the distribution Vi^f) of 7 = log g /(g). Ev- 
idently, log7-'(7) = 7 + logP(exp7). Therefore, an 
assumption P(g = 0) = c ^ 0, implies V = 7 + logc 
for 7 — > —00. 

Fig 1. shows clearly that log ^(7) decreases more 
quickly than 7 for all ensembles we consider. This guar- 
antees that P(g) — > as g — > 0. Let us note that it is 
almost impossible to obtain last result form the studies 
of P(g) on the linear scale [(Tof . 

The small- g behavior of P(g) is easy to estimate also 
from the distribution P(z\) of the smallest parameter 
Z\ . Indeed, small values of g require large values of Z\ . 
Neglecting contributions of other channels, we have 

1 f 2e 1 f°° 

- ^ P(g)dg = - J_ P{z 1 )dz 1 (2) 

with s = exp— ii. In the limit e — > the integral on 
the LHS reads ~ P(g), g = s. RHS could be found an- 
alytically for special form of P(zi). In particular, for 
Wigner surmises P{z\) — n/2{zi) 2 x zicxp(— 7r/4 x 
[zi/(zi)] 2 ) we obtain that P(g) - ^-i-constxlogg / 2 
with const = 47^2 ■ Consequently, P(g = Q) = 0. 
Fig. 2. assures that P{z\) decreases more quickly than 




Figure 2: Probability distribution of (normalized) z\ 
for 4D and 3D Anderson model. Solid line is Wigner 
surmises Pw{z) — § zexp[— j z 2 ]. For the mean value 
(zi) see Table 1. 

Wigner surmise for large z\ in orthogonal ensemble for 
both 3D and 4D systems. This assures that P(g) — > 
as g — > 0. 

Linear behavior of the distribution P{z{) for small 
z\ (see left side of Fig. 2) guarantees nonzero prob- 
ability that the first channel is fully open. Indeed, 
if P{z\) ~ C x z\ for Z\ — > 0, then the probability, 
that the first channel contribution to the conductance, 
gi = 1/ cosh 2 (zi), equals to 1, is C. This explains the 
origin of the characteristic bump in the distribution 
P{g) for g = 1. In Fig. 1, the bump is clearly visible 
for both 3D and 4D systems. 

Fig. 1 (right) confirms that P(g) decreases more 
quickly than exponentially for large g. This is easy to 
understand on the basis of the analysis of the statis- 
tics of z's presented in ||. Fig 3. shows mean values 
and variances of some smallest z's for both 3D and 4D 
system. Evidently, (zi) ~ 0(1) and variances var Zi 
decreases quickly with index i. Consequently, the con- 
tribution to the conductance from the second (higher) 
channel is, due to ([!]), small (negligible). To estimate 
this contribution, we note that all higher Zi, i > 2, are 
normally distributed. |^]. Their mean and variances 
has been estimated as (z^ ~ (zi) xi l ^ d ~ 1 ^ and var Zi ~ 
{z i )~( d ~ 2 ' ) [|l6). Although this result holds only in the 
quasi-one dimensional limit, where the mutual correla- 
tions of z's are negligible, they serve as a good quanti- 
tative estimation also for true d-dimensional cubes. As 
i is seen in Figure 3, this agreement is better for 4D 



Figure 3: Var zj as a function of (zi) (only for (z) < 15) 
for 3D (open symbols) and 4D (full symbols) orthog- 
onal systems. Note the system-size invariance of pre- 
sented parameters (at least for i < L). 

than for 3D,. Then, the probability to find g ~ n is 
less than exp[— (z n )/2varz„] ~ exp[— const x n d ^ d ^ 1 ^ > ] 
and 

P{g) ~ exp -const x c^" 1 ). (3) 

We conclude that presented numerical data for 3D 
and 4D Anderson model prove the system size invari- 
ance of the conductance distribution at the critical 
point. Although the distribution depends on the di- 
mension and symmetry of the system, we found its 
common features, namely exponential decrease of P(g) 
for g > 1, and a decrease of P(g) to zero for g = 0. We 
show that the form of P(g) can be analyzed on the 
basis of the statistics of parameters z introduced by 
relation (1). This analysis is more simple for higher 
dimension, where the statistical correlations of zs are 
supposed to be less important. 
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